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Abstract 

It was recently observed in a lattice QCD measurement that the chiral con- 
densate in the quenched approximation shows dramatically different behavior 
in the three Zs-equivalent de-confined phases. We argue that this phenomenon 
can be understood qualitatively as an effect of Z3 twists on fermionic fields. 
Quarks under these Zs-twists become global anyons and hence display differ- 
ent thermodynamic properties. We further show that the lattice data can be 
roughly modeled by a Nambu-Jona-Lasinio type Lagrangian with a minimal 
coupling to a constant gauge field ^0 = 27rnT/3 (with n = 0, ±1), which 
arises naturally from the non-trivial phase of the Polyakov line. 
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I. INTRODUCTION 



Understanding chiral symmetry breaking and its restoration at finite temperatures in 
QCD is of great interest. Many lattice simulations have been performed to study the ex- 
pected phase transitions [|l|. However, the underlying dynamics responsible for these phe- 
nomena are still unclear. Hence, it is sometimes useful to explore regions of QCD parame- 
ters beyond their physical domain. One such effort in lattice QCD is to study the fermionic 
Greens functions in the background dynamics of pure gauge fields. In such a study one 
ignores all the effects of internal quark loops. At low temperatures no essential difference is 
observed between such studies and those that take into account the effects of quark loops. 
However, at finite temperatures these effects can change the physics qualitatively. In the 
pure gauge theory there is a first order phase transition between a confined and a de-confined 
phase signaled by a non-vanishing expectation value of the Polyakov line. This is related to 
a global Z3 symmetry breaking On the other hand, inclusion of massless quarks breaks 
this global Z3 symmetry explicitly and chiral symmetry restoration becomes relevant, with 
the chiral condensate as the order parameter. The precise nature of the chiral transition 
depends on the number of fiavors appearing through internal quark loops. 

Interestingly, in lattice QCD one can measure a quenched chiral condensate even after 
ignoring internal quark loops from the dynamics. Such a condensate has been measured 
in the past near the de-confining transition. A first order transition was observed in the 
quenched chiral condensate at the de-confining temperature [Q. However, no efforts were 
made to distinguish between the behavior of quarks in the three de-confined phases labeled 
by an element of Z3. Recently, in an effort to understand the effects of internal quark loops 
near the chiral phase transition, this issue was revived by a careful and precise measurement 
of the quenched chiral condensate 0. It was observed that, above the de-confining transition, 
the chiral condensate crucially depends on the Z3 phase in which the gauge dynamics settles. 
Even the chiral symmetry restoration appears to occur at different temperatures, depending 
on the phase. This somewhat surprising result needs an explanation. In this article we 
present a qualitative understanding of the physics underlying these observations. 

In a somewhat different context the effect of the boundary condition in the Gross-Neveu 
model, with one of the dimension compactified, was studied recently |^ . It was found there 
that the fermion, which obeys the usual antiperiodic boundary condition, becomes a global 
anyon when a non-trivial twist is introduced via a constant background gauge field. The 
anyon interpretation is supported by an explicit decomposition of the fermion propagator 
into a sum over winding numbers in the compactified direction. As a consequence of this 
statistics transmutation the thermodynamic properties of the Gross-Neveu model become 
sensitive to the value of the twist. We suggest that a similar physical setting is realized in 
quenched QCD in the various de-confined phases. The non-trivial phase of the Polyakov line 
provides a constant background gauge field in the temporal direction and turns quenched 
quarks into global anyons. 

In section 2, we show that in QCD one can define a chiral condensate with a Z3 twist in 
general. The usual chiral condensate is the one where the twist is zero. The other two are 
new quenched order parameters of chiral symmetry breaking. We then argue that the chiral 
condensates measured in are related to these Z3 twisted condensates in the quenched limit. 
In section 3, we elucidate the meaning of the Z3 twisted chiral condensate by considering an 
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effective action which governs these quenched observables. In section 4, we show that the 
quahtative features of the lattice data can be obtained by the Nambu-Jona-Lasinio model 
in the large Nc limit, minimally coupled to the suggested constant gauge field. In section 5, 
we present some observations and conclusions. 



II. Z3 TWISTED CHIRAL CONDENSATES 



The chiral condensate in QCD is the order parameter for chiral symmetry breaking and 
is defined formally as 

(?^) = ^ [ [dA] tr—^ Det{D + m,) exp(-5,[A]) (1) 

Zj J U -\- Tfiq 

where D[A) = 7^((9^ — iA^^) is the Dirac operator and rriq is the quark mass and Sg is the 
gauge action. Z is the partition function. The above definition can be made more precise 
on the lattice. 

We now introduce a chiral condensate with a Z3 twist 6, given by 

(CO. = ^ f [dA] tr—^ DetiD + ruq) exp(-5,[A]) (2) 

where 



D{A)e = ^^{d^ - iA^)+ieT^^ (3) 

with 6 = 0, ±^ and T the temperature. The mass in the observable, rriQ, is deliberately 
chosen to be distinct from the usual quark mass, rriq. Thus {(() is a quenched observable 
and can be measured on the lattice even in the presence of dynamical fermions. The physical 
content of this new observable will be elucidated shortly. Notice that Dg is simply the original 
Dirac operator with the temporal gauge field shifted by a phase of Z^. Further note that 
(CC)e = (CC)-e- We can show this by noting that D{A) and D{—A) have the same spectra 
due to a charge conjugation symmetry. However, in general we will have (CO 221 7^ (COo- 
This is related to the fact that fermions break the Z3 symmetry. 

In general, the dynamical quarks in eq. (^, will always pick the phase in which the 
gauge field distribution, {^4^}, makes the expectation value of the Polyakov line P real and 
positive, i.e. 

(P(x)) = (^^Tr I^P JKm^AoM dt)^"^ = \{P)\ (4) 

The gauge field distribution of a Z3 rotated phase would be {A^} + 9 T S^^. Now the 
Polyakov line is given by \{P)\ e*^. The relevant Dirac operator for measuring the chiral 
condensate is D{A)g defined in eq. ^. Thus measuring the usual chiral condensate in the 
various Z3 phases is naturally equivalent to measuring (CC)e the quenched limit of the 
de-confined phase. 

The effect of similar 6 terms on the thermodynamics has been studied in the past for 
simple local fermionic field theories 0. There it is argued that the 6 term acts like a twist 
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in the boundary condition for the fermions and converts them into global anyons. This 
conversion then affects the critical temperature dramatically. We suggest that the same 
physics also happens in the present context. We anticipate that, if we study {C()e as a 
function of the temperature, it is likely that the condensate with ^ = and 6 = ±^ will 
behave differently. In particular if there is a phase transition in the rric^ —>■ limit it is quite 
likely that it will occur at different temperatures depending on 6. 

The above discussion then suggests that it is likely that at some temperature in the de- 
confined phase (CC)o = whereas (CC)^ O5 as observed in 0. However, in the confined 
phase the gauge dynamics is Z3 symmetric in the quenched limit. The Z3 twist in this phase 
ceases to be meaningful, since the gauge field ensemble is no longer clustered around any of 
the three 9 values there. It is of course still possible to study the Z3 chiral condensates, as 
defined in eq. 0, even in the presence a dynamical quark mass rriq, because of the relative 
shift in the Dirac operator in the observable and in the internal quark determinant. In the 
presence of finite rriq the Z3 symmetry is broken and thus it is possible to track the three 
phases into the confined phase and study their evolution towards the chiral limit. We do 
not know if such a study is physically relevant since the Z^ condensates are in general some 
quenched order parameters. In the remainder of this article we will be interested only in the 
quenched limit. 



III. AN EFFECTIVE ACTION FOR THE TWISTED CONDENSATES 

In this section we suggest a local field theory which governs the behavior of the Z3 twisted 
chiral condensates. Formally one can rewrite eq. (^, introducing fermionic {() and bosonic 
(0) ghost fields in the spirit of 0, 

(COe = \ I CC exp \-S,[A] + ^ (D + m,)V^ + C{De + m^)C + 4>\De + m^)^] (5) 

This enlarged theory contains another chiral symmetry in the limit — ^ for arbitrary real 
quark mass rriq. It is now obvious that the order parameter associated with this symmetry 
is (CO- 0^16 can then formally integrate over the quark fields ip, the ghost fields (p and the 
gauge fields to generate an effective action for the ( fields. Thus 

(CO. = exp [-S^Zo,m„T] (C, C) + Cb'd, + T 7° + m^)C] (6) 

We have allowed for the possibility that the effective action could depend on the temper- 
ature through the couplings apart from the usual anti-periodic boundary conditions in the 
Euclidean time direction for the ( fields. 

In general S^^l e m^T](CX) would be arbitrarily complicated and possibly contains non- 
local interactions. It seems hopeless to understand the full structure of the effective action. 
However, it may be interesting to obtain some qualitative insight into the essential physics 
involved by studying a simplified model that is motivated phenomenologically and at the 
same time is tractable. 

In the quenched limit, the chiral symmetry between fermionic Greens functions is gov- 
erned only by the dynamics of the gauge fields. Non-perturbative gauge fields can in general 
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produce chiral symmetry breaking. One immediate consequence is that the quenched chi- 
ral condensate need not vanish. This feature, in principle, need not have any connection 
with confinement. Thus it is not unnatural to assume that, for some temperature range, 
chiral symmetry breaking happens even in the deconfined phase. Consequently we are led 
to study chiral symmetry breaking in the absence of confinement. To this end, we assume 
that the effective interaction generated after integrating out all other fields is dominated by 
the lowest dimensional term, 

SHUmd'^^ = ^ [(C • C)^ - (C ■ 750^] (7) 

with A being some cutoff scale specified later. Thus in this model the 9 dependence enters 
only through the quadratic term shown in eq. ^. We expect this term to produce the 
difference between the complex and the real phases as in [Q. 

As mentioned earlier, in the confined phase there is no difference between the various 
twisted condensates. This is due to the Z3 symmetry of the dynamics. Thus the effective 
action must have no knowledge of 9. This is in accordance with the fact that the expectation 
value of the Polyakov line is zero in the confined phase, it is meaningless to include a 9 
term in the confined phase. Here no attempt is made to construct an effective field theory 
that describes the confinement physics also in a natural way. However, as long as we are 
interested only in a model for the chiral condensate, it may be sufficient to drop the 9 
term from the effective action in the confined phase. Additional effects of gauge dynamics 
are mimicked by a temperature dependence of the coupling G{T). Further since the de- 
confining transition is a first order transition, an abrupt jump is possible in G{T) at the 
de-confinement temperature. 



IV. MODELING OF LATTICE DATA 



To have a semi-quantitative understanding of the lattice data we start with the La- 
grangian 



^(C,C) = c- 

or equivalently 

^(C, = (■ 



GiT) 



(C-C)'-(C-75C)^ 



2G{T) 



(8) 



(9) 



Since this model is non-renormalizable we consider a momentum cutoff in the spatial 
momentum |p| = A. We will solve the model in the large limit and fix the parameters 
in such a way that the model reproduces the lattice data obtained in [^] qualitatively. The 
effective potential in the large limit is given by 



NcA^ 
2G{T) 



{2ny 



(10) 



where 



4 



Un = [{2n - 1)71 + 9]T 



(11) 



The Matsubara sum can be carried out by the standard contour integral technique, yielding 

G{T) , . , \. 1 1 



G{T) 



a 



7r2A2 



(a + m) / dp — 

Jo -C/„ 



1 - 



^Ep/T+iO _|_ ]_ ^Ep/T-ie _|_ I 



12) 



with Ep = + (cr + m)2. In the above equation the anyon behavior is clearly displayed 
from the modified thermo-distribution, which interpolates smoothly from the Fermi-Dirac 
{9 = 0) to Bose-Einstein {6 = it) distributions. 

The gap equation in the large Nc limit reduces to 



0" 



G(T) 



cr + m M p dp 
7r2A2 io "eT 



1 - 



1 



1 



^Ep/T+iB _|_ I ^Ep/T-ie _|_ I 



(13) 



The full effective potential can be obtained after integrating eq. [12 



Kfr(^) - Kfr(O) 



_ G'(T)T 
G(r) I y " 7r2A2 



/ p^ dp In 
Jo 



cosh(wp2 -|- (cr -)- rriY/T) + cos 6* 



cosh(-\/p2 _|_ rn? /T) + cos 6' 



(14) 



We refer the reader to |Q for a detailed discussion of the effects of 9. Here we only try to 
model the lattice data of 0. 

The parameters in the NJL model are fixed in a conventional way [||. Assuming the 
cutoff A to be about 1 GeV, we set A/T^ = 5.5 and G{T) /{2ti'^) = 1.163 for T < in 
order to reproduce the chiral condensate at low temperatures with Tc around 200 MeV. The 
9 term is dropped in the confined phase for reasons already explained. The current quark 
mass on the lattice {rrii^), and in the effective action (m) can in general be different. Here we 
assume^ m = Z ttiq. Then Z = 1.48 can be cleanly determined by fitting the condensate 
{9 = 0) in the chirally symmetric phase. However, in order to reproduce the lattice data 
quantitatively a non-trivial T-dependence in G(T) needs to be introduced in the deconfined 
phase. The first order nature of the de-confining transition is incorporated by a jump in 
G(T). Then a simple quadratic form G(T) /{2tt)^ = 1.385 — 0.385(T/Tc)^ is sufficient for 

T>Te. 



In addition to contributions proportional to m^, m could get other contributions, such as from 
the Ua{1) anomalies. We will ignore such additive contributions here. 
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FIGURES 
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FIG. 1. This graph shows the lattice data. The behavior of the condensate in the large 
Nc Nambu-Jona-Lasinio model described in the text, produces the solid lines. The lattice data 
corresponds to a = 0.001, where a = with Nt = 4. 

The lattice data corresponds to A^i = 4 and the lattice spacing is given by a = ^rj^- 
Then the value of the condensate computed from the NJL model should be compared with 
the lattice data as 

-W^^-^^ = = gaKtJ = J-G^'^ ^^^^ 

where a is obtained from the solution of the gap equation. 

The data from lattice simulations and the model are given in figures |l| and ^ [| 
It is seen there that the NJL model, with its parameters properly tuned, is capable of 
describing the data semi- quantitatively. However, we would like to stress the fact that the 
final outcome of the NJL model sensitively depends on the precise temperature dependence. 
A small change in the numerical value of the coupling G can change the evolution of the 
complex chiral condensate drastically. Thus in this regard we find the model somewhat 
unnatural, even though it is possible to capture the qualitative physics. 

In figure ^ the current mass dependence of the condensate is compared between the 
model and the lattice data at T/Tc = 1.067. At this temperature the complex phase lattice 



^ The actual data was given in terms of the lattice coupling. This has been converted to dimen- 
sionless less units using the critical temperature as the scale. The bare lattice coupling /? = 5.692 
was used as the critical coupling for the de-confinement temperature. 
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data begins to show interesting power law dependence, and the real phase lattice data is 
essentially linear 0. The real phase data at T = 1.067 was used to fix the constant Z and 
hence will agree well with the model. However, it is obvious that the present model will not 
reproduce the power law behavior of the real phase lattice data at lower temperatures [^. 
Similarly we do not expect the model to reproduce the power law behavior of the complex 
phase either. In fact we do not know if any universality arguments apply to these data. 
However, we do see the qualitative difference between the real and the complex phases. 
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FIG. 2. This graph compares the behavior of the condensate as a function of the current 
quark mass at T = 1.067Tc in the large Nambu-Jona-Lasinio model described in the text with 
the lattice data. This is approximately the temperature where a possible transition occurs in the 
complex phase, (r) refers to = and (c) refers to = ib27r/3. 



V. DISCUSSION AND CONCLUSIONS 

In this work we have clarified the precise meaning of the chiral order parameter considered 
in a recent lattice study. The observed Z3 phase dependence of this order parameter in 
quenched QCD in the deconfined phase has been qualitatively explained. The relevant 
physical picture is that the non-trivial Z3 twists in the temporal boundary condition, due to 
the survival of expectation value of the Polyakov line, turn quarks into global anyons. The 
qualitative behavior of the lattice data can be modeled simply by a temperature dependent 
NJL model with a minimal coupling to a constant gauge field = QT with ^ = and 
±27r/3. 

For m = 0, the NJL model would predict a continuous phase transition when Q = 
±27r/3, in addition to the jump associated with the first order de-confining phase transi- 
tion. The present lattice data are consistent with this prediction. The second transition 
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temperature is estimated to be about 5 to 10 percent higher than the de-confining transition 
temperature. 

As this work was being completed we received an article |T^, which discusses the same 
issue along very similar lines using the NJL model. However, based on our calculation, we 
would like to point out that, on a quantitative level the model does not refiect one important 
feature of the gauge dynamics. It is clear from figure |l| that the lattice data for the complex 
phase evolves rapidly above Tc towards a possible second transition at about T = 1.07(3) Tc. 
In the above model such a rapid evolution required very fine tuning of the coupling G{T). 
Thus it is very likely that the gauge dynamics is not naturally captured in this model. In fact 
the rapid evolution just mentioned is also seen in other physical quantities like the entropy 
density [ill] same region. Thus an obvious direction to extend the present work, is to 

obtain a better model including the gauge dynamics. This may help in understanding the 
reason for the rapid evolution seen just above T^. 

As was suggested in section ||, the Z3 condensates as defined by eq. ^ can be studied even 
in the presence of small dynamical quark masses rrig. As quenched order parameters these 
may be measured in lattice simulations. However, more theoretical work would be necessary 
to see if they are interesting. Finally, in this article we have concentrated on the lattice 
data for SU{3) gauge theories. A similar discussion in the case of SU{2) leads to interesting 
possibilities since the de-confining transition is second order in this case. Extensions to 
SU {N) with > 3 may also be interesting and has been attempted in [ffO . 
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